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Abstract 

We employ 3D Langevin Dynamics simulations to study the dynamics of polymer chains translo- 
cating through a nanopore in presence of asymmetric solvent conditions. Initially a large fraction 



I (> 50%) of the chain is placed at the cis side in a good solvent while the trans segments are placed 

in a bad solvent that causes the chain to collapse and promotes translocation from the cis to the 
trans side. In particular, we study the ratcheting effect of a globule formed at the trans side created 



' by the translocated segment, and how this ratchet drives the system towards faster translocation. 



Unlike in the case of unbiased or externally forced translocation where the mean first passage time 
(r) is often characterized by algebraic scaling as a function of the chain length N with a single 
! scaling exponent a, and the histogram of the mean first passage time P{t/ (r)) exhibits scaling, we 



find that scaling is not well obeyed. For relatively long chains we find (r) ^ N°' where a ~ 1 for 
e/ksT > 1. In this limit, we also find that translocation proceeds with a nearly constant velocity 
of the individual beads(monomers), which is attributed to the coiling of the globule. We provide 
an approximate theory assuming rotational motion restricted on a 2D disc to demonstrate that 
there is a crossover from diffusive behavior of the center of mass for short chains to a single file 
translocation for long chains, where the average translocation time scales linearly with the chain 
length A*". 

PACS numbers: 87.15.A-, 87.15.H-, 36.20.-r 
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INTRODUCTION 



Attempts to understand the dynamics of viral invasion and infection dynamics of 
DNA and other biopolymers passing through porous media {i.e. cell membranes), and 
direct medical applications such as gene therapy and drug delivery, have rendered the field 
of polymer translocation a very active field of research in recent years . Much work has 
been done to understand the physics involved in the translocation process. Analytical work 



by Sung and Park Muthukumar |^|, Chuang, Kantor and Kardar [7|, l8|, Dubbeldam and 
coworkersfl. Pauja a„d cowoAersQ. a„d others^, supplemented by a vast amount 
of numerical work [1^-|21|, has brought profound physical insight into the problem at hand. 
In particular, much has been learned about making translocation faster in a controllable 
fashion, as this should be beneficial in biological systems. 



To this end, a process known as Brownian ratcheting 22|] was discovered early on and 
discussed further in Refs. 23|]- 27|]. Brownian ratchets are mechanisms by which translocation 
can be driven, or carried out more efficiently 27|] . As discussed by Simon, Peskin, and 



Oster 



221 ]. thermal or chemical asymmetries in the system can be used to extract useful work 



(i. e. translocation of the polymer) from the thermal bath in accordance with the Second Law 

n 

of Thermodynamics [27|. As the polymer translocates, it experiences considerable back and 
forth motion due to thermal fluctuations. If the part of the chain that is on the trans side is 
modified in such a way as to prevent backward motion through the pore, its random motion 



22|. 



will be biased and translocation through the pore is notably faster [24|. This modification 
to the chain which causes a biased translocation is often called a Brownian ratchet 
A Brownian ratchet can manifest itself in many different ways. There can exist bindin. 
particles that bind as chaperones on the trans side 
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28|, glycosylation can be used |24j . 



or the chain can be ti ght ly bound into a coil on the trans side via some method, usually by 
having a bad solvent 16[ or reducing the solvent's pH 24 • ^ bad solvent, the polymer 
chain undergoes a coil-globule transition to form a highly interacting spherical-like polymeric 
configuration with a radius of gyration Rg that scales with the number of monomers N as 



(Rg) ~ A^i/^y |3o|. 



In the present work, our aim is to study the infiuence of ratcheting on the dynamics of 
polymer translocation as induced by an asymmetry in solvent quality between the cis and 
trans compartments. We focus on the case of two-sided translocation [l^, where a fraction 
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of the polymer chain is initially placed at the trans compartment, with good solvent on the 
cis side and bad solvent on the trans side. This asymmetry in the solvent condition induces 
a bias in the entropic barrier controlling translocation in such a way as to effectively drive 
the polymer to the trans side. To characterize the dynamics of translocation, we analyze the 
waiting times, velocities, and effective forces on the individual monomers inside the pore. 
We find that for long chains with high attraction strengths, e/kBT, waiting times vary only 
slightly until the last beads emerge, at which point entropic effects become dominant. Since 
the velocity of the beads is inversely proportional to the waiting time, we get roughly constant 
velocity in this regime. This is apparent only for long chains with relatively high interaction 
strengths; for shorter chains with lower interaction strengths, the center of mass velocity of 
the polymer introduces an dependence into the velocity that cannot be overlooked when 
calculating how the average translocation time scales with A^. We used this idea to make an 
approximate estimate for the dependence of the average translocation time. We find for 
short chains the translocation exponent a ~ 2, while in the large limit, we find a — ?■ 1. 
We end our analysis by analyzing the histograms of the mean first passage time(MFPT). We 
note that since a varies between 2 and 1, we do not see universal scaling for a given value of 
e/kBT. However, in the large A^ limit, we begin to see scaling manifest itself clearly, with a 
scaling exponent approaching the predicted value of unity. 



II. MODEL 



We use the Langevin equation to study the Brownian motion of particles in solution. It 
is a statistical, stochastic differential equation of the form for each bead i: 

mri{t) = -WUi - Vv,{t) + Wi(t), (1) 

where the total interaction, 



Ui = U 



FEME 



LJ ' 



(2) 



is the sum of the finitely extensible nonlinear elastic(FENE) spring potential interaction 31 | 
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UpENEirij) = --kRlln 



Rq 



and the Lennard- Jones interaction between neighboring monomers. 



ULj{rij) = iSi 



a. 
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cr. 



(3) 



(4) 



The term W(t) describes the influence of Markovian white noise due to the solvent, which 
is not taken into account explicitly here. It satisfies the fluctuation-dissipation relation 

(W(t) ■ W(r)) = QkBTT6,j6{t - r). (5) 

To model asymmetric solvent conditions on the cis and trans sides, interaction cut-off values 
were modified using a cut-off matrix. Each particle was given a label (either 1 or 2, depending 
on whether it was on the cis side or trans side, respectively). The cutoff values for the trans 
side were set higher (rc,22 = 2.5cr) than those for the cis side {r^^n = rc,i2 = i^c,2i. 



2^/^(t). This ensures that the monomers on the trans side interact with a "bad" solvent 321 ] . 
while the monomers on the cis side are in a good solvent characterized by a self avoiding 
random walk with Flory exponent u = 0.588 {Rg ~ iV^) in 3D. The different cut-off values 
introduce a chemical potential difference, Afi, between the compartments. We will study 
the corresponding solvent asymmetry for various monomer-monomer interaction coupling 
strengths, e. In particular, we will study how this solvent quality asymmetry drives the 
system towards a much faster translocation. 

The purely repulsive wall consists of one monolayer of immobile L J particles of diameter 
1.5cr on a triangular lattice at the xy plane at 2; = 0. The pore is created by removing 
the particle at the center. The reduced units of length, time and temperature are chosen 
to be cr, cry^, and e/ks respectively. For the spring potential we have chosen k = 30 
and Rij = 1.5a, the friction coefficient F = 1.0, and the temperature is kept at 1.5/kB 
throughout the simulation. 

For a chosen fraction of the monomers at the cis /trans we equilibrate the chain for a time 
on the order of the Rouse relaxation time r ~ A^^"*"^*^, where the Flory exponent u = 0.588 
in 3D . The chain is then allowed to translocate using a time step of dt = 0.005. As the last 
bead exits the pore, a translocation event is completed and the process repeated for 2000 
times for averaging. 

One way to understand the dynamics of a polymer chain translocating under such highly 
asymmetric conditions, is to study the analytic form of its free energy. Following Muthuku- 
mar, the free energy for m translocated monomers is given by 5| 

^ = (1 - 7^ ln(m) + (1 - 70 HN - m) + m^. (6) 

Here, 7'=0.5, 0.69, and 1 for Gaussian, self-avoiding, and rod-like chains, respectively. Driv- 
ing force is easily obtained from the free energy by differentiating with respect to monomer 



index: 



ksT dm m {m — N) ksT 

This is the driving force of the system in units of fc^T for particular values of 7J, 72, and 
A/i. As we shall demonstrate below, there is a delicate balance between the frictional and 
driving forces that will tend to set the system at a constant velocity. 

The translocation time for a chain is a function of the number of monomers on the trans 
side at the beginning of the translocation process Ntr{t = 0). If we start our simulation 
having 50% of the chain on the trans side, corresponding to the "two-sided" translocation 
first considered in Ref. [l^, i.e. Ntr{t = 0)/N = 0.5, this corresponds to releasing the chain 
down a downhill entropic barrier and therefore, the probability for successful translocation, 
P{Ntr{t = 0), should be unity, which is indeed the case in our simulation. This probability 
decreases drastically as the fraction Ntr{t = 0)/N is less than 0.5, and especially for long 
chains, the probability for a successful translocation is very small as shown in Fig. 1 for a 
chain of length = 64. In the present work, we have studied the cases for Ntr(t = 0) = 0.5A^ 
and Ntr{t = 0) = 0.25A^. By comparing these two sets of data, we have extracted the results 
for the limit Ntr{t = 0) — 0. As expected, we recover uniform scaling of the probability 
distribution for the MFPT with a translocation scaling exponent a — )■ 1 in the large A^ limit. 

III. RESULTS AND DISCUSSION 

To get an idea of the translocation process, we show typical snapshots of a translocating 
chain in Fig. 2(a) and Fig. 2(b) at different stages of the translocation process. At t = 0, 
the fraction of the chain that is located on the cis side is characterized by the equilibrium 
Flory exponent u ^ 0.588, that corresponds to the good solvent condition. Since the trans 
)art of the chain is in a poor solvent and the temperature is below the 6-temperature 



32|, it will form a globule which is expected to grow as a function of time. Comparing 
the snapshots for e/fefiT = 0.5 (Fig. 2(a)) and e/ZcsT = 1.5 (Fig. 2(b)), we note that 
the globule formed by the translocated segments becomes more compact as the strength of 
the interaction increases. We have checked the A^ dependence of the radius of gyration for 
chains immediately after the translocation process as shown in Fig. 3. For larger interaction 
strength, we find Rg ~ A^'^'^'^, which is consistent with the A^ dependence of a compact 
spherical globule. It is worth mentioning that if all the globules for different chain lengths 
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were perfect spheres and fully relaxed, then Rg ~ ]\f0.33_ p^j, ejkBT = 1.5, we note from the 
snapshots that the spheres are very compact and hence a dependence of Rg ~ A^'^-^'' implies 
that the spheres formed by the translocated segments for different chain lengths are close 
to equilibrium. For e/ksT = 0.5, the corresponding exponent extracted from the slope of 
= 64, 128 and 256 yields Rg ~ A^'^'^^, which is less than 1/3. This indicates that the trans 
side of the chain does not have sufficient time to relax during the translocation process as 



discussed earlier in the literature 



11 



19 



We have monitored several quantities during the translocation process. First, we have 
analyzed the waiting time as a function of monomer index (Fig. 4), where we have defined 
waiting time to be the total time each bead spends at the pore divided by the total translo- 
cation time (r) for the whole chain to cross the pore {i.e. W{m) = {rim)) / (r)), where r(m) 
is the total time bead m spends at the pore). The notation (■■) indicates ensemble average 
over 2000 iterations. We notice that W{m) initially increases and then decreases to almost 
zero, and the position of the maximum increases with the chain length A^. We also find that 
the peak position is an increasing function of the interaction strength e/ksT. Previously, 
waiting time of a monomer for a translocating chain was studied in great detail |l4|. For a 
homopolymer undergoing externally forced translocation, the residence time increases and 
becomes maximum for a monomer index rrimax > N/2; it then decreases (more rapidly than 
the rise) almost linearly, the position of the maximum being skewed towards m > N/2. 
When an attractive interaction is present for the translocated segments, the barrier that the 
monomers at the cis side are pulled through is effectively skewed. However, this is different 
from applying a force only on the monomer inside the pore. We note that data in Fig. 4 
are similar to the case of forced translocation. However, for the last ~10% of monomers, 
the residence time decreases very rapidly. This becomes more pronounced with increasing 
strength of the attractive interaction. It is evident from the residence time plots that longer 
chains with larger attraction strengths result in residence times that vary only slightly over 
the trajectory of the translocation, until the last beads emerge. We note that the waiting 
time of each monomer is inversely proportional to the velocity of that particular monomer 
at the pore. We have plotted the reciprocal of the waiting time function (scaled by the 
appropriate factors) for a chain length of A^ = 128 (Fig. 5). We note that the velocity and 
inverse waiting time function collapse onto the same graph and are relatively constant up 
to the last few beads, at which point they both drastically increase due to entropic forces. 
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We will explore the velocity at the pore for all the chains lengths below. 

The behavior of the residence time when translated to the average velocity of the monomer 
at the pore for various chain lengths = 16 — 256 and for e/ksT = 0.5, 1.0 and 1.5 is shown 
in Figs. 6 (a)-(c). As predicted from the residence time plots, the velocity for long chains 
is virtually constant for the entire translocation process until the end, when the last beads 
emerge from the pore. We note that the driving force should be proportional to the velocity. 
From Eqn. (7) given for the driving force, we note that the term involving the reciprocal of 
the difference N — m, blows up as the last few beads emerge from the pore, explaining why 
velocity increases for the last few beads at the pore. This is most evident for longer chains 
with higher e/kBT (Fig. 6(b)-(c)). 

To further understand our results, we have studied the force experienced at the pore 
as a function of the monomer index Fig. 7(a)- (c). In order to get a better idea about the 
interaction of the chain with different solvents on either side of the pore we have not shown 
the force arising out of the high frequency phonons from the anharmonic spring potential 
and have shown only the LJ contribution to the force in our plots. For e/ksT = 0.5 — 1.5 and 
relatively long chains, we see a rather flat force curve, which is close to F = 0, in agreement 
with V = const discussed earlier (Fig. 7(a)). We interpret this result as being indicative of 
a force balance between frictional and driving forces Ff and F^r, respectively, at the pore. 
The driving force F^r is given by in Eqn. (7.) For force balance to occur, we must then 
have the following condition: 

{F,r) ^ (Ff). (8) 

In the large limit, the driving force will be governed primarily by the chemical potential 
difference between the compartments A/j,. In this limit 

{Fdr{N ^ oo)) ^ Aj2^T{V{N ^ oo)). (9) 

Thus, we see that the coiling velocity Vc will be proportional to the chemical potential 
difference: 

(fc) ~ ~ const. (10) 

The main result is that in the large A^ limit, the driving force and velocity are independent 
of A^ and only dependent on the chemical potential difference between the compartments. 
This driving force is exactly balanced by the friction experienced by the monomers, which is 
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proportional to the velocity of the beads. Thus we see that in the large N limit the condition 
^ F ^ and v ^ const is well obeyed. 

In the case of short chains, the force becomes negative gradually for the last beads, 
whereas for long chains the change is rather drastic (Fig. 7(a)- (c)). We interpret this negative 
force as being a result of the last few beads, (having a relatively large velocity) while still 
being at the cis compartment but in the vicinity of the pore, getting absorbed by the globule 
on the trans compartment, resulting in a deceleration as they escape the cis compartment. 
This results in a negative force on the last few monomers. 

Next we present a simple scaling argument for estimating how the average translocation 
time should scale as a function of A^. It is based on the observation that the translocation 
dynamics corresponds to a "coiling" of the chain around the collapsed globule on the trans 
side with coiling velocity Vc (Fig. 8). If we assume that the attractive force on the 
bad solvent side is directed towards the center of the globule, we can write Fg ~ mv^/R, 
which gives Vc ~ (RFe/m)^^'^ . For the collapsed globule close to equilibrium R ~ N^^^ 

1 /2 

and m N ^ and thus v^. ~ (^N~'^/^F^ . Due to spherical symmetry, we expect that the 
force will be proportional to the number of monomers in a disk of radius R ~ N^^'^. Thus, 
Fe ~ A^^/^, from which we extract that Vc = const. Now, {t)v ~ i?^, where R^ ^ N and 
V = 1/N + Vc = 1/N + const. The dependence comes from the velocity of the center of 
mass, and will only contribute in the low limit. Thus we can see that in the low limit, 
the 1/A^ term dominates and thus (r) ~ A^^, which is supported by our results below. On 
the other hand, for large A^, the constant will dominate, thus arriving at the large A^ limit 
of (r) ~ A^, also supported below. Thus for single file translocation induced by coiling, we 



would expect a scaling exponent close to a ~ 1, much lower than 3D 
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breed translocation 



2l|. 



where a ~ 1.37 — 1.6 depending on the rate of translocation p/7|. 

We have tested the scaling argument above for several values of e by studying the his- 
tograms of the MFPT for several values of chain length, namely A^ =16, 32, 64, 128, and 
256 to see how they scale as a function of A^. Following previous work 20| we have used 
nonlinear regressions of the form f{x) = Ax^ exp{—Cx). The maxima for these curves occur 
ai X = B/C. We have used the position of the maxima for each chain length to obtain the 
mean first passage time (r) from which we can extract the scaling exponents given in Table 
I. Scaled MFPT histograms are shown in Fig. 9(a)-(c) for e/kbT = 0.5 — 1.5. We notice from 
Table I that the scaling exponent a decreases when extracted from successive larger values 
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of A^. If we look at the spectrum of exponents calculated for short and long chains and for 
various values of the interaction strength e qualitatively we notice that for short chains and 
weaker e the exponent reflects diffusive behavior, while for larger combination of and e 



the exponent is 
chain lengths 



17 



ess than the corresponding exponent for forced translocation for similar 



20 



21| . In the next section we provide theoretical argument why the 



translocation behavior from short and low e is dominated by diffusive behavior while 
for long and large e, this diffusive behavior crosses over to a "single file" translocation 
asymptotically reaching a scaling exponent a — 1 in this limit. Evidently, for this reason we 
do not see data collapse of the scaled histogram for the MFPT accros the board. However, 
it is worth noticing that for long A^ and large e {N = 128 and 256 and e/kBT = 1.5) we 
notice a near perfect data collapse. One can see this trend from Eqn. (7). For short chains 
we see that chain length plays quite an important role in determining the driving force on 
the polymer. For relatively long chains, as we showed earlier, this A^ dependence is washed 
away, and the only contribution to the driving force is the chemical potential difference. 
Thus for N ^ oo the driving force becomes independent of chain length. 

We have also compared how the MFPT data collapse on a single master curve when 
we use the initial condition, Nt{t = 0) = N/4 instead of Nt{t = 0) = A^/2 as shown in 
Fig. 10. We notice that data collapse and scaling is more closely obeyed in this regime. 
The scaling exponent for e/ksT = 1.5 continues to suggest a tendency towards a — )■ 1.0 
in the large A^ limit, as is also predicted by Wei et al. in their studies of the effect of 



solvent quality asymmetries on the translocation process [16|, although studied differently 
and using a different model for the solvent conditions. Their studies indicate that polymers 
translocating under different solvent qualities have a scaling law that varies from (r) ~ N'^^'^^ 
to (r) ~ A^, which is quite close to our results. Our present results are also consistent 
with Muthukumar's analytical expression for the translocation time as a function of A^ for 
various conditions for the chemical potential difference j^. Calculations in Muthukumar's 
work show that for symmetric barriers, the translocation time will scale as (r) ~ A^^. For 
asymmetric barriers and long chains, if the entropic terms in the free energy equation are 
small compared to the term involving A/i, then the translocation time scales linearly as 
(r) ~ A^ for A^||A/i|| > 1, and scales as (r) ~ A^^ for A^||A/i|| < 1, which is consistent with 
our present results. For relatively short chains [i.e. N = 16) and weak coupling strength 
e, we note a scaling exponent close to a ~ 2. For longer chains with stronger coupling 
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constants e, we notice a trend towards unity, in complete agreement with Muthukumar's 
analytic expression. We expect that for longer chains, this limit will be even more closely 
reached due to the single file nature of the translocation induced by the coiling of the 
globule on the bad solvent compartment. This is also reflected in in Fig. ll(a)-(c) we show 
the translocation time (r) plotted against on a log-log scale. In each attempt to scale the 
MFPT histograms, we have used the linear fit slope of the entire log-log plot (^.e (r) ~ A^"). 
The exponents support our claim that the the translocation interpolates from diffusive to a 
single file behavior. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated using 3D Langevin Dynamics simulations the properties of a ho- 
mopolymer translocating through a nano-pore with a solvent asymmetry. In our model, 
there is good solvent on the cis side of the pore, while on the trans side there is bad sol- 
vent. This creates an effective driving force on the polymer and leads to the emergence of 
a collapsed globule on the trans side during the translocation process. We have used a free 
energy argument to show that the driving force is relatively insensitive to chain length in 
the large N limit and is governed mainly by the chemical potential difference between the 
compartments. As is evident from our force plots, there is a delicate balance between this 
constant (in the large A^ limit) driving force and the friction force experienced by the beads. 
Consistent with this idea, we find that the velocity of the beads at the pore is relatively 
constant in the large A^ limit, which is attributed the constant coiling velocity occurring on 
the "trans" compartment. Furthermore, we note that scaling is not well obeyed in the low A^ 
limit as is evident from our mean first passage time histograms. We interpret this as being 
a consequence of the A^ dependence of the driving force (and thus velocity!) in the low A^ 
limit. For longer chains, however, we note that the driving force(and velocity) become quite 
insensitive to changes in the chain length, and we retrieve scaling with a scaling law close to 
our predicted large A^ limit of (r) ~ A^. This can be interpreted as a crossover phenomenon 
from diffusive type translocation to a single-file driven translocation. Our studies might be 
relevant for translocation of biopolymers accros cell membranes. 
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FIGURE CAPTIONS 

Fig. 1: Probability of successful translocation as a function of the initial number of 
monomers on the trans side for a chain of length = 64. 

Fig. 2(a): Snapshots of a translocating chain of length = 256 for e/ksT = 0.5 at times 
(a) t = 0, (b) 0.25r, (c) 0.5r, (d) 0.75r, and (e) l.Or respectively. 

Fig. 2(b): Snapshots of a translocating chain of length A^ = 256 for e/kBT = 1.5 at times 
(a) t = 0, (b) 0.25r, (c) 0.5r, (d) 0.75r, and (e) l.Or respectively. 

Fig. 3: Variation of (Rg) as a function of number of translocated monomer Ntr (log-log 
plot) for e/kBT = 1.5 (black circles) and e/kBT = 0.5 (red squares) respectively. 

Fig. 4: Average waiting time as a function of monomer index normalized by the maximum 
waiting time for various chain lengths. The symbols circles (black), squares (red), diamonds 
(green), triangle up (blue), and triangle left (magenta) correspond to the chain lengths 
A^ = 16, 32, 64, 128, and 256 respectively (color online) for (a) e/kBT = 0.5. (b) 
e/kBT = 1.0, and (c) e/kBT = 1.5. 

Fig. 5: Inverse of waiting time (red diamonds) and velocity of the monomer beads 
(blue circles) plotted as a function of monomer index for chain length A^ = 128 and 
e/kBT = 0.5. We note that the two graphs almost fall on top of each other and are rela- 
tively constant until the end of the translocation, when the last beads emerge out of the pore. 

Fig. 6: Average velocity on the monomer beads inside the pore as a function of monomer 
index for chain lengths A^ = 16, 32, 64, 128, and 256 respectively for (a) e/kBT = 0.5, (b) 
e/kBT = 1.0, and (c) e/kBT = 1.5. The symbols have the same meaning as in Fig. 4. 
We note that in the large A^ limit, the velocity of the monomers becomes almost constant 
excepting for the last few monomers. 

Fig. 7: Average force on the monomer beads inside the pore as a function of monomer 
index for chain lengths A^ = 16, 32, 64, 128, and 256 respectively for (a) e/kBT = 0.5, (b) 
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e/ksT = 1.0, and (c) e/ksT = 1.5. The symbols have the same meaning as in Fig. 4. To 
analyze the force, we have taken only the LJ contribution (neglected the smallest the fluctu- 
ations). In agreement with the velocity plots, we notice that the force at the pore is not only 
constant, but also zero in the large limit, which agrees with the constant velocity that we 
find at the pore, indicating a force balance at the pore between friction and driving force. 
We also notice that the force drastically becomes negative as the last few beads translocates. 

Fig. 8: Schematic of a polymer coiling ideally around a collapsed globule. The force is 
directed towards the center and is proportional the number of monomers in a 2D disc 
around which the translocating polymer coil. 

Fig. 9: Scaled histograms for the MFPT for different chain lengths A^=16, 32, 64, 128, and 
256 for (a) e/Aj^jT = 0.5, (b) e/kBT = 1.0, and (c) e/ksT = 1.5 for Ntr{t = 0)/N = 0.5. 
The symbols have the same meaning as in Fig. 4. Here, we have used the linear fit slope for 
the entire Log-Log plots (Fig. 11) for e/Zc^T = 0.5 — 1.5 as scaling exponents. In general for 
these values of e/ksT we do not notice universal scaling for all chain length N. However, 
for longer chains and larger attraction strengths (inset in (b) and (c)), we begin to see 
scaling emerges quite clearly for a ^ 1. 

Fig. 10: (a) Unsealed and scaled MFPT Histograms for chain length = 64 (green 
diamonds) and = 128 (blue triangle- ups) with Ntr/N{t = 0) = 0.25 as the initial 
condition, (a) for s/ksT = 0.5 where we notice a scaling exponent close to the value 
attained for the corresponding case of Ntr/N(t = 0) = 0.50 shown in Fig. 9. (b) For 
e/ksT = 1.5; we notice a scaling exponent closer to the theorized large limit of unity. 
We also note that scaling is more closely obeyed in this case. 

Fig. 11: Variation of (r) as a function of A^ (log-log plot) for (a) e/Zc^T = 0.5, (b) e/ksT = 
1.0, and (c) e/ksT = 1.5 respectively. The local values of the slope (a) are indicated in the 
graph. 
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Table I: Effective scaling exponents ((r) ~ N") for various combinations of and 
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